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Abstract. The NA60 experiment has measured muon pair production in In-In collisions at 158 AGeV 
at the CERN SPS. This paper presents a high statistics measurement of meson production. 

Differential spectra, yields, mass and width are measured as a function of centrality and compared to 
. , ■ previous measurements in other colliding systems at the same energy. The width of the rapidity distribution 

' is found to be constant as a function of centrality, compatible with previous results. The decay muon polar 

angle distribution is measured in several reference frames. No evidence of polarization is found as a function 
of transverse momentum and centrality. The analysis of the pr spectra shows that the (j) has a small radial 
flow, implying a weak coupling to the medium. The Te// parameter measured in In-In collisions suggests 
that the high value observed in Pb-Pb in the kaon channel is difficult to reconcile with radial flow alone. 
The absolute yield is compared to results in Pb-Pb collisions: though significantly smaller than measured 
by NA50 in the muon channel, it is found to exceed the NA49 and CERES data in the kaon channel at 
any centrality. The mass and width are found to be compatible with the PDG values at any centrality and 
at any pr- no evidence for in- medium modifications is observed. 

PACS. 25.75.Nq Quark deconfinement, quark-gluon plasma production, and phase transitions - 25.75.-q 
Relativistic heavy-ion collisions - 25.75.Dw Particle and resonance production - 14.40.Cs Other mesons 
with S = C = 0, mass < 2.5 GeV - 12.38.Mh Quark-gluon plasma 



1 Introduction 

QCD predicts the occurrence of a phase transition from 
hadronic matter to a deconfined plasma of quarks and glu- 
ons, when sufficiently high energy densities are reached. 
Strangeness enhancement was proposed long ago as a sig- 
nature of this phase transition [Tl[5] . The meson is a key 
experimental probe in this context: due to its ss valence 



quark content it permits to measure strangeness produc- 
tion. It was suggested that, being a hidden strangeness 
state, the (j) should not be sensitive to canonical suppres- 
sion, and can thus be used to test strangeness enhance- 
ment models [1I1[51[M7]. 

Moreover, it has been argued that the spectral func- 
tion of the (j) or of its decay products could be modified 
in the medium. This would reflect in changes of its mass 
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and partial decay widths. Mass shifts of vector mesons, 
as a consequence of the change of the quark condensate 
in the medium, were originally proposed by Brown and 
Rho [8]. Subsequently, the meson mass dependence on 
the medium density was calculated with QCD sum rules 0. 
Modifications of the width in the context of this dropping 
mass scenario were investigated, for instance, in the model 
of Ref. \lOi . As an alternative approach, both the mass and 
width changes were studied in the framework of hadronic 
many body models (see [H] for an up-to-date discussion) , 
which predict a broadening and no or small mass shift. 
Finally, it was suggested that in-medium changes could 
also appear in cold nuclear matter [T2lll3j. 

It is important to point out that, contrary to the p 
meson which has a lifetime of 1.3 fm, much smaller than 
the lifetime of the fireball (7-10 fm), the has a lifetime of 
44 fm. Therefore in-medium effects on the mass or width 
might be difhcult to isolate experimentally, and might only 
appear at very low px- A possible way of tracing down 
those modifications is the comparison of the particle yields 
and spectra in the dilepton and in the kaon channel: since 
the (j) is almost at the threshold for the cj) KK chan- 
nel, even a small change in its spectral properties should 
lead to a significant difference in the two channels [14]. A 
discrepancy in the two channels could also be ascribed to 
final state effects related to in medium kaon absorption or 
rescattering [l5lfT6]. 

The KEK-PS E325 experiment reported evidence for 
a large excess to the low mass side of (/) peak in p-Cu at 
12 GeV for (3^ < 1.25, corresponding to slow ip mesons, 
having a higher chance to decay inside the Cu nucleus 
(the effect was not observed in p-C and for higher values 
of P"f) [17j . This excess was explained in terms of an in- 
medium modified component with mass decreased by ^ 
35 MeV and width broadened by a factor 3.6 with respect 
to the PDG value. The reported yield associated to this 
excess is very significant: 20% of the total (f) yield. 

(f> production was studied extensively at the SPS by 
several experiments. NA49 measured (j) production in p-p, 
C-C, Si-Si and Pb-Pb at 158 AGeV in the K+R- 
channel [TSlfTM^ . These data show that the (j) yield breaks 
Npart scaling and that the enhancement, quantified as 
{(/)) /Npart, saturates already for Si-Si central collisions. 
Pb-Pb collisions were also studied in different experimen- 
tal conditions by the NA50 experiment, which found con- 
trasting results. NA50 studied the nfi channel f2T] 
with acceptance limited to pr > 1-1 GeV, while the NA49 
measurement in the </> K^K~ channel was limited to 
Pt < 1.6 GeV because of statistics. Both experiments ob- 
served an enhancement of the (p yield with the size of the 
collision system. The absolute yields, however, disagree 
by a factor ~ 4 in full phase space [52]. A recent NA50 
analysis of a new sample confirmed previous findings [23] . 
The NA49 and NA50 results disagree also on the T^f f pa- 
rameters extracted from the transverse mass spectra fitted 
with the thermal ansatz. The NA49 T^ff values are con- 
sistently larger and show a stronger centrality dependence 
than those measured by NA50. While this could be partly 
due to the radial expansion, which flattens the distribu- 



tions at low rriT, the effect is too large to be explained by 
radial flow alone. 

Naively, kaon absorption or rescattering could prevent 
the reconstruction of in-matter —^ KK decays, in par- 
ticular at low transverse momentum, while the (j) mesons 
decaying in the lepton channel would not be affected. This 
would lead to reduced yields and enhanced T^ff in the 
KK channel as compared to the fip. channel. Meson mod- 
ification in the medium could also lead to differences in 
the two channels, as mentioned above. The difference be- 
tween NA49 and NA50, however, is larger than predicted 
by models including all these effects piIllTilfBllTBtei] . The 
discrepancies in absolute yields and Tg// became known 
as the (f) puzzle. 

Production of in central Pb-Pb collisions was re- 
cently studied also by CERES [5S] both in the kaon and 
dielectron channels. The measured T^f f and yield are com- 
patible with the NA49 results. However, due to the large 
uncertainty in the dielectron channel, the authors could 
not rule out completely a possible difference in the two 
channels, though significantly smaller than the NA49/NA50 
difference: the maximum possible enhancement in the lep- 
tonic over hadronic channel was estimated to be 1.6 at 95% 
confidence level. 

In this paper we report on new measurements in the 
muon channel done by the NA60 collaboration in In-In col- 
lisions at 158 AGeV incident beam energy, which can help 
to clarify (j) production in nuclear collisions. Transverse 
momentum and yields have been studied as a function of 
centrality. Transverse momentum spectra, measured in the 
range — 3 GeV, and yields are compared to previous re- 
sults, providing further insight on the puzzle. The mass 
and width of the are studied as a function of centrality 
and Pt, in order to search for possible in-medium effects. 
The results are compared to previous measurements and 
existing theoretical predictions for In-In collisions. In ad- 
dition, the rapidity and the decay angular distributions in 
three reference frames (Collins-Soper, Gottfried- Jackson 
and the helicity frame) are also measured. In elementary 
collisions, the (j) decay angular distributions can convey 
information on the production mechanism [2&. No mea- 
surements exist in p-p, while the AGGMOR collaboration 
performed measurements of the polarization in pion, 
proton and kaon - Be interactions [27ll28] . A small but 
significant transverse polarization was reported at low pT 
in the Gottfired- Jackson frame, which was related to the 
parton fusion model to estimate the average mass of the 
fusing partons. More recently, NA60 measured the angu- 
lar distributions in In-In collisions, integrated in central- 
ity and for pT > 0.6 GeV, for the excess dimuons below 
1 GeV (dominated by tt+tt" ^ p — > /i+/i~) and the vec- 
tor mesons uj and |29j . The reported absence of any 
anisotropy indicates the completely random orientation 
of the interacting hadrons in 3 dimensions as expected in 
the case of a thermalized medium. In this paper we ex- 
tend these results by performing the study as a function 
of centrality and px . 
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Fig. 1. The NA60 experimental apparatus. 
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Fig. 2. Detail of the vertex detectors layout. 



2 The NA60 apparatus 

In this section the detectors of the NA60 apparatus used 
in the analysis are briefly described. A more detailed de- 
scription can be found in Ref. [30] . 

The 17-ni-long muon spectrometer, previously used by 
the NA38 and NA50 experiments [21] , is complemented by 
new silicon detectors in the vertex region. The experimen- 
tal apparatus is shown in Fig. [1] The muon spectrometer 
is composed of 8 multi-wire proportional chambers and 4 
scintillator trigger hodoscopes, divided into a "forward" 
and a "backward" telescope by the 5-m-long ACM (Air 
Core Magnet) toroidal magnet. A 120 cm thick iron wall, 
located just before the last trigger hodoscope and after 
the tracking chambers, ensures that only muons can trig- 
ger the apparatus, without degradation of the tracking 
accuracy. The angular acceptance is 35 < < 120 mrad; 
the rapidity coverage for the </> is 3 < y < 4.2. 

A 12 \i hadron absorber, mostly made of graphite, is 
placed in front of the tracking system of the muon spec- 
trometer, as close as possible to the target to minimize 
the fraction of hadrons that decay into muons. The iden- 
tification provided by this "muon filter" leads to a very 
selective trigger, enabling the experiment to run at very 



high luminosities. However, the material in the hadron 
absorber also degrades the kinematics of the muons, be- 
cause of energy loss fluctuations and multiple scattering. 
This affects in particular the angular resolution, which is 
'-^ 10 mrad for typical muons from decays. 

The layout of the vertex region is shown schematically 
in Fig.[2l Downstream of the target and before the hadron 
absorber, a silicon tracking telescope [32] embedded inside 
a 2.5 T magnetic dipole field measures the charged tracks 
produced in the interaction. The muons in the vertex re- 
gion are then identified by matching them to the tracks 
reconstructed in the muon spectrometer. As described in 
the next section, this significantly improves the perfor- 
mance of the detector. The telescope used during the in- 
dium run consisted of 16 independent detector planes, pro- 
viding 11 tracking points, arranged along the beam axis 
over a length of ~ 26 cm, starting at ~ 7 cm from the tar- 
get centre. The telescope covers the angular acceptance of 
the muon spectrometer, with the first planes smaller than 
the last ones. Each plane consists of several single-chip 
pixel detector assemblies, 96 in total, which are mounted 
on a planar ceramic support. The assemblies are made 
of radiation tolerant, 750 /im thick ALICElLHCb pixel 
readout chips [33J, bump-bonded to 300 [ivn thick p-on-n 
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Fig. 3. Longitudinal coordinate of interaction point. 

silicon pixel sensors, with 425 x 50 fim^ pixels. Each of the 
single chip assemblies is a matrix of 32 columns and 256 
rows. They have been shown to remain fully functional up 
to a radiation dose of 12 Mrad [SJ- The detector coped 
well with the radiation level of '^l Mrad per week, reached 
in the cells around the beam axis of the first planes in the 
NA60 environment [321 . The telescope tracks all charged 
particles produced in the collision, thus providing an esti- 
mate of the collision centrality. 

Incoming ions are tracked by a silicon microstrip beam 
tracker, composed of two stations placed 20 cm apart 
(30 cm and 10 cm upstream of the target). The 400 /im 
thick sensors have 24 strips of 50 fim pitch, allowing to 
infer the transverse coordinates of the interaction point, 
at the target, with a resolution of 20 fim 36 . They are 
operated at cryogenic temperatures (130 K), to remain 
sensitive to the beam ions even after collecting consider- 
able radiation doses [37]. The "back planes" of the four 
sensors extend over a larger area and constitute a very 
good beam intensity counter. 

Embedded inside the hadron absorber, a Zero Degree 
Calorimeter (ZDC) provides an alternative estimate of the 
collision centrality by measuring the energy carried by the 
projectile nucleons which have not taken part in the in- 
teraction (spectators) [38ll39| . The ZDC and the beam 
tracker were used to measure the beam intensity in the 
absolute yield measurement (sec. 14. 4p . 

3 Data sample and analysis procedure 

During the 5- week- long run in 2003, 230 million muon pair 
(dimuon) triggers were recorded, with two different set- 
tings of the ACM magnet current (4000 A and 6500 A). 
The results presented in this paper refer to the 4000 A 
data sample, which amounts to ~ 50% of the total statis- 
tics, as the 6500 A setting suppresses the low mass accep- 
tance and it is optimal for the study of the J/tp. 

The data were collected with a 158 AGeV indium 
beam, incident on 7 indium targets of 17% total inter- 



action probability. The targets were placed in vacuum in 
order not to pollute the data sample with In-air interac- 
tions. The z coordinate of the interaction vertex is shown 
in Fig.[3l The resolution on the position of the interaction 
vertex is 10 — 20 /im in the transverse coordinates and 
~ 200 /im in the longitudinal coordinate. The target iden- 
tification is extremely good, with negligible background 
between different targets. Other elements present along 
the beam line are also clearly visible in the picture, like 
the target box windows and the beam tracker stations. 

In the following subsections we illustrate briefly the 
most relevant aspects of the reconstruction and selection 
of events, background subraction and acceptance correc- 
tion. 



3.1 Event reconstruction and selection 

As a first step, the muon tracks are reconstructed in the 
muon spectrometer. About half of the triggers are dis- 
carded either because it was not possible to reconstruct 
two muons or because at least one of the muons traversed 
the 4-m-long iron poles of the ACM. These tracks suffer 
from a very strong multiple scattering, which would sev- 
erly degrade the momentum resolution. The muon tracks 
reconstructed in the muon spectrometer are then extrapo- 
lated back to the target region and matched to the tracks 
reconstructed in the vertex spectrometer. This is done 
comparing both their angles and momenta, requiring a 
matching Isss than 3. Once identified, the muons are 
refitted using the joint information of the muon and the 
vertex spectrometers. These tracks will be referred to as 
matched muons. Muon pairs of opposite charge are then 
selected. The matching technique drastically improves the 
signal-to-backround ratio and the dimuon mass resolution, 
which decreases from ~ 80 MeV (using only the informa- 
tion from the muon spectrometer) to ^ 20 MeV at the 
(j) peak, independent of centrality [40] . 

Events with more than one interaction vertex or with a 
(reconstructed) charged track multiplicity smaller than 4 
are discarded |40j . The matched dimuon mass spectrum is 
shown in the left panel of Fig. [U The same picture shows 
also the two background sources: combinatorial and fake 
matches. These are discussed in detail in the next subsec- 
tion. As it can be seen, the dominant background is the 
combinatorial background - the fake match background is 
almost one order of magnitude smaller. The w and (j) peaks 
are very well resolved and the rj peak is also visible. The 
right panel of Fig. |3] shows an expanded view of the uj and 
<j> region after the subtraction of backgrounds. The sample 
consists of ^--^70 000 4> events. The signal/background ra- 
tio, integrated over centrality, is 1/3 below the (f> peak, 
ranging from about 3 in peripheral events to about 1/5 in 
central ones. 

To study (j) production, events in a small mass window 
centred at the 4> pole mass (0.98 < M < 1.06 GeV) are se- 
lected from the mass distribution after the combinatorial 
background subtraction. The continuum below the (f> is 
accounted subtracting the events in two side mass win- 
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Fig. 4. Left panel: raw mass distribution (continuous line), combinatorial background (dashed line), fake tracks (dot-dashed 
line) and mass distribution after combinatorial background and fake tracks subtraction (continuous curve with error bars) . Right 
panel: mass region around the 4> peak and mass windows used in the analysis after backgrounds subtraction. The dot-dashed 
distribution is the total fake matches distribution. 



dows (0.88 < M < 0.92 GeV and 1.12 < M < 1.16 GeV), 
shown in the right panel of Fig. [H 

No specific cuts in transverse momentum, rapidity or 
CoUins-Soper angle are applied. As already mentioned, the 
detector introduces an intrinsic rapidity cut 3 < y < 4.2 
at the mass. 



3.2 Background subtraction 

The dimuon sample is affected by two sources of back- 
ground: the combinatorial background and the fake mat- 
ches. The former is the contribution of uncorrelated muon 
pairs, mostly coming from the decay of pions and kaons, 
and is subtracted with an event mixing technique PUIHT] . 
The NA60 apparatus triggers also on like sign pairs IJ,~ 
and ^"'"/j,''"), which, at SPS energies, contain only uncorre- 
lated muons. Two muons from different like-sign events are 
randomly paired to build a dimuon sample which is uncor- 
related by construction. The comparison of the real and 
mixed spectra for the like-sign pairs leads to an estimate 
of the accuracy of the background subtraction of ~ 1% 
over the full dimuon mass spectrum. The like sign sample 
is also used to determine the absolute normalization of 
the combinatorial background [101I¥T| . The combinatorial 
background is shown as a dashed line in the left panel of 
Fig. 12 

The background from fake track matches comes from 
the fact that, when the multiplicity is high, a muon track 
can be associated to more than one track in the vertex 
telescope with an acceptable matching x^, and the track 
with the lowest may not be the correct match. It can 



also happen that the correct track in the vertex telescope 
in absent altogether. The contribution of the fake matches 
can be estimated in two different ways. The first approach 
is an overlay Monte Carlo technique, where a Monte Carlo 
dimuon is reconstructed on top of a real event. Selecting 
the match with the smallest matching it is possible to 
determine the fraction of fake matches as a function of 
mass, transverse momentum, etc. The left panel of fig. [S] 
shows the matching for the correctly matched tracks 
and for the fake matches, as obtained with an overlay 
Monte Carlo simulation. The stricter the matching cut, 
the higher the signal to background ratio will be. However, 
a tight cut reduces significantly the statistics. The second 
method is based on the event mixing technique, which ex- 
tracts the probability distributions of fake matches from 
real data alone. The basic idea is to match the tracks in 
the muon spectrometer from one event to the vertex tracks 
of a different event [40ll41j . All the matches obtained in 
this way are fake by construction. The two methods agree 
within 5%. 

The right panel of Fig.[5]shows the Monte Carlo dimuon 
invariant mass for cj) events when one or both muons are 
not correctly matched, for signal events and for all recon- 
structed pairs. This gives an idea how incorrect matching 
degrades the kinematics, resulting in a broadened (j) peak 
with a ~ 100 MeV. Notice however, that this is not the 
total fake matches background under the peak, because 
also other processes can contribute - continuum from open 
charm or tails from the p—cu region. The total fake matches 
distribution is shown in detail as a dashed line in the right 
panel of Fig. [H This shows that the fake matches, in the 
specific case of the (f> meson, can also be directly elimi- 
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Fig. 5. Matching (left) and mass distribution (right) of signal and fake matches for simulated cf) events. 



nated by the subtraction of the side windows, without the 
need to estimate them with the methods outlined above. 

In the following, the analysis is repeated subtracting 
the fakes both with the event mixing technique and with 
the side windows. 



3.3 Centrality selection 

All results are studied as a function of the collision cen- 
trality. Events are divided into five centrality bins by us- 
ing the charged tracks multiplicity as summarised in Ta- 
ble [T] The reconstructed minimum bias multiplicity dis- 
tribution depends on the target because of acceptance, 
reconstruction efficiency and secondary particles contami- 
nation. The distributions in the seven targets are rescaled 
so that they overlap to that in the first target. The latter 
is then fitted with the multiplicity distribution obtained 
with a Glauber Monte Carlo simulation, assuming that 
the number of reconstructed tracks is proportional to the 
number of participants, with a smearing due to the reso- 
lution and inefhciencies of the detector. It is assumed that 
the effects due to acceptance, reconstruction efficiency and 
secondaries are embedded in the proportionality constant 
and smearing factor. The distributions and mean values 
of the number of participants Npart and of the number 
of binary collisions Ncou , are obtained applying the same 
multiplicity cuts used in the real data to the simulated 
data. Alternatively, the real data can be corrected for ac- 
ceptance, reconstruction efficiency and secondaries using a 
Monte Carlo simulation [42,43]. The number of measured 
tracks can be thus converted to the number of primary 
particles event by event, as a function of target and cen- 
trality. The distribution of primaries is then matched to 



Table 1. Centrality bins used in the analysis, selected using 
the charged tracks multiplicity. The intervals used in the five 
centrality bins are reported in the table, together with the 
mean values of the number of participants (A^part) and of the 
number of binary collisions (Acoii) (see text for details). 
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the distribution of Npart obtained by a Glauber Simula- 



tion. We find N, 



part 



dN,h/dv\s,r MM- 



The difference on the mean Npart and N^oii values for 
each bin obtained with the 2 methods is taken as an es- 
timate of the systematic uncertainty and is found to be 
about 10% for the peripheral bins and 5% for the central 
bins. 



3.4 Acceptance and reconstruction efficiency 
correction 

After subtraction of the combinatorial background, the 
physical continuum and the fake tracks, the data are cor- 
rected for acceptance and reconstruction efficiency. These 
are estimated with an overlay Monte Carlo simulation, us- 
ing the real hit occupancy in the detector superimposed 
to Monte Carlo muon pairs. The same cuts used for the 
real data were applied to the simulated data. An overall 
95% efficiency was assumed for the pixel detector. The in- 
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as a function of pr for the (j). 



as a function of y for the 4> for < pr < 0.5, 0.5 < pr <1 and 
Pt > 1 GeV. 



put sources' kinematics are generated using the Genesis 
event generator [3S] , tracked through the apparatus using 
GEANT 3.21 and reconstructed using the NA60 analysis 
framework. 

The study of the differential spectra would require a 
multi dimensional acceptance correction in rapidity, trans- 
verse momentum and Collins-Soper angle. Since the low 
populated bins in the phase space corners can introduce 
large errors, a fiducial cut should be applied, at the cost 
of reducing the effective statistics. To avoid this loss, a 
one dimensional correction is applied to each kinematic 
variable under study, integrating over the remaining ones. 
This requires a careful tuning of the Monte Carlo input 
kinematics which was adjusted with an iterative procedure 
to the measured distributions for each of the five central- 
ity bins. To validate this procedure, the full analysis chain 
(event selection, reconstruction, acceptance correction, fit- 
ting) is applied to Monte Carlo events. It is verified that 
the corrected kinematic distributions and the estimated 
parameters are consistent with those used as input in the 
Monte Carlo. 

The acceptance times reconstruction efficiency for the 
as a function of px is shown in Fig. [H] The apparatus 
has a good acceptance down to zero pT, due to the addi- 
tional dipole field in the vertex region. Indeed, low mass 
and low pT dimuons, which would not be accepted by the 
muon spectrometer, are bent into the acceptance by the 
vertex dipole magnet, as opposed to NA50, which used 
the same muon spectrometer, but had no acceptance for 
Pt < 1.1 GeV. Further details on the rapidity and an- 
gular distribution acceptance will be given in the specific 
subsections on the results. 



3.5 Systematic effects 

The main sources of systematic error are the uncertainties 
on the input parameters of the Monte Carlo simulation 
(which reflect on an uncertainty in the acceptance) and 
on the estimate of the physical background below the 0. 
The former is estimated by changing the values of the 
parameters in the input Monte Carlo within \a of the 
measured value. The latter, by varying the side windows 
offset (from 50 MeV to 70 MeV) and width (from 30 MeV 
to 50 MeV) . The contributions from these two factors are 
added in quadrature. Uncertainties in the combinatorial 
background and in the fake matches subtraction play a 
small role here, as the peak is very well resolved. 

Several consistency checks were performed to test for 
the presence of hidden systematic effects [46] : the match- 
ing cut was varied between 2 and 3, the kinematic cuts 
2.9 < ?/ < 3.9 and |cos0cs| < 0.5 were applied, the analy- 
sis was repeated independently in the 4 ACM/PT7 fields' 
combinations, a 2-dimensional acceptance correction (vs. 
y andpr) was used and the different methods for fake sub- 
traction were implemented. Most of those checks did not 
lead to significantly different results, with the exception 
of the matching cut variation. This leads to a small 
but not negligible difference which is added in quadra- 
ture to the systematic error. This difference is most likely 
due to the uncertainty in the energy loss compensation 
or to a small residual misalignment in the pixel telescope 
planes. The results on the yield, mass and width are af- 
fected by further sources of systematics, due to the pa- 
rameters needed for their computation. Details are given 
in sec. 14.41 and in sec. 14.61 The number of participants es- 
timated with the Glauber model is affected by a 10% un- 
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Fig. 8. Left panel: Gaussian a of the rapidity distribution of the (p meson, as a function of centrality; full circle: (p fijj, in In-In 
(NA60); full upward triangle: (p K+ K' in p-p (NA49); full downward triangle: K+ K' in C-C and Si-Si (NA49); open 
circle: (j) K~ in Pb-Pb (NA49). Right Panel: Rapidity distribution of the cj) meson; open symbols are obtained reflecting 

the measured points around mid-rapidity. 



certainty in the peripheral bins, which decreases to ~ 5% 
in the most central, as explained in sec. 13.31 



4 Results 
4.1 Rapidity 

The rapidity acceptance for < pt < 0.5, 0.5 < pt < 1 
and Pt > 1 GeV is shown in Fig. [71 The apparatus has 
acceptance in the interval 3-4.2 depending on the target. 
There is a small correlation y — pr such that low y tends 
to select larger px- The rapidity is studied as a function 
of centrality and the width estimated with a gaussian fit 
in the range 3.2 < y < 4.1. The midrapidity value is fixed 
to 2.91. The width as a function of centrality is shown in 
Fig. [8] (left panel) and within the errors one can deduce 
either a constant or a slightly decreasing trend. The rapid- 
ity distribution integrated in centrality is shown in Fig. [8] 
(right panel). The resulting width is 1.13 ± 0.06 ± 0.09. 
This is in agreement with previous NA49 measurements 
in other collision systems at the same energy ^19u20j . as 
shown in Table [21 



4.2 Angular distributions 

In general, the angular distribution of the decay muons in 
the (f) meson rest frame can be related to the elements of 
the spin density matrix. 



Table 2. Gaussian width of the rapidity distribution of the 
(j) meson in several collision systems, at 158 AGeV incident 
beam energy. The second quoted error for the NA60 point is 
the systematic error. 



Experiment 


System 


ay 


NA49 


p-p 


0.89 ± 0.06 


NA49 


C-C 


1.16 ±0.10 


NA49 


Si-Si 


1.27 ±0.10 


NA60 


In-In 


1.13 ±0.06 ±0.09 


NA49 


Pb-Pb 


1.22 ±0.16 



In this paper the inclusive dN/d cos 9 distributions are 
studied as a function of transverse momentum and cen- 
trality. The polar angle 9 is the angle between the /i+ and 
the z axis as measured in the following three widely used 
reference frames. In the helicity frame, the axis is defined 
as the direction of flight of the (j) meson in the target- 
projectile rest frame. In the Gottfried- Jackson frame (GJ), 
the axis is defined as the projectile direction in the rest 
frame \47_^. Finally, in the CoUins-Soper frame (CS) the 
axis is given by the bisector of the projectile and the neg- 
ative target direction in the rest frame [48] . 

In all the three frames, the (j) acceptance as a function 
of cos 9 has a strong dependence on transverse momen- 
tum. In particular, the apparatus has zero acceptance for 
high I cos 1 and low pT- The acceptance vs cos 9 for the 
Pt windows < pt < 0.5, 0.5 < pt < 1 and px > 1 GeV, 
is shown in Fig. [9l for the CS and GJ frames. The most 
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Fig. 9. Geometrical acceptance times reconstruction efficiency as a function of cos 9 for the for < pt < 0.5, 0.5 < pr < 1 
and Pt > 1 GeV in the CS frame (left) and in the GJ frame (right). 



critical situation occurs for the CoUins-Soper angle, while 
the helicity frame is similar to the GJ frame. 

The dN / d cos distribution is fitted with the function 
1 + A-cos^ 6*, where the parameter A is related to the degree 
of polarization. The spectra, integrated in centrality and 
PT: are shown in the left panel of Fig.lTOl The polarizations 
obtained from the fits are zero within errors in all cases: 
Xhel = -0.07 ± 0.06 ± 0.06, Agj = -0.07 ± 0.06 ± 0.06, 
Acs = -0.13 ±0.07 ±0.06. 

The polarization is then studied differentially in the 
three px windows < pt < 0.5, 0.5 < pt < 1 and > 1 
GeV, without centrality selection. No significant devia- 
tions from the pt integrated values are found. 

The results as a function of centrality and integrated in 
Pt are shown in the right panels of Fig. 1101 The polariza- 
tion is compatible with zero in the helicity and Gottfried- 
Jackson frames, independent of centrality. In the CoUins- 
Soper frame there is a hint for a negative polarization in 
the most peripheral bin. However, this is within 2a com- 
patible with zero and then consistent with the absence of 
any effect as seen in the other frames. 

These results as a function of pt and centrality are con- 
sistent with the results integrated in centrality with px > 
0.6 GeV reported in a previous analysis by NA60 [551. 
They show the absence of polarization, supporting the fact 
that particles are produced from a thermalized medium. 
The absence of polarization extends down to peripheral 
events, suggesting that it could be zero also in elementary 
collisions at this energy. 



4.3 Transverse momentum 

The acceptance-and-efficiency corrected transverse mass 
distributions are shown in the left panel of Fig. [11] for 
the different collision centralities, normalized in absolute 
terms. 

The transverse momenta arc first fitted in the range 
< < 3 GeV with the thermal ansatz: 

I dN I dN , , , 

= — -1— exp(-mT/Te//). (1) 

rriT drriT Pt dpT 

The effective temperature Tg/ / as a function of the num- 
ber of participants is depicted in the right panel Fig. 1111 
In general one can notice that the Te// shows an initial 
rise followed by a saturation in the most central bins. 

The exponential form assumes thermal emission from 
a static source. In heavy ion collisions, this is not correct 
due to the presence of radial flow, but eq. ([T]) was never- 
theless used by several experiments to fit the data. Since 
the rriT distributions are not purely exponential. The Te/ / 
extracted from an exponential fit will depend in general 
on the fitted px range. It can be shown that, at first order, 
the Teff value at low px is given by T^/f '-^ Tq + to/3|,, 
where Tq is the temperature at freeze-out, m is the parti- 
cle mass and (3x is the average transverse velocity of the 
expanding source [49] . At high px, on the other hand, the 
Teff extracted from the exponential fit is approximated 
by T - Tov/(1 + /3t)/(1-/3t) ■49|. Thus, what is mea- 
sured is a Px dependent effective temperature, larger than 
the freeze-out temperature because of the radial flow. 

The influence of radial flow is investigated by repeat- 
ing the flts in two sub-ranges < pT < 1.6 GeV (low 
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Fig. 11. Left panel: (p transverse mass distributions in indium-indium collisions as a function of centrality; from top to bottom: 
central to peripheral spectra. Right Panel: Centrality dependence of the T^ff parameter; fit performed over the full NA60 range 
(0 <PT < 3 GeV). 
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Table 3. Ccntrality dependence of the Te// parameter in diflterent fit ranges. The second quoted error is the systematic error. 



(A^part) 



15 
41 
78 
133 
177 



T,ff (MeV) 
< pT < 3 GeV 



T,ff (MeV) 
< pT < 1.6 GeV 



T,ff (MeV) 
1.1 <pt < 3 GeV 



209 ± 4 ± 2 
232 ± 3 ± 3 
245 ± 2 ± 4 
250 ± 2 ± 3 
249 ± 3 ± 4 



206 ± 5 ± 2 
229 ± 5 ± 2 
253 ± 5 ± 9 
258 ± 5 ± 5 
253 ± 8 ± 11 



24 

78 
133 
177 



228 ± 6 ± 4 

237 ± 5 ± 1 
245 ± 4 ± 5 
244 ± 5 ± 3 



Pt) and 1.1 < pr < 3 GeV (high pr), which also corre- 
spond to the experimental windows covered by the NA49 
and NA50 experiments, respectively. Since the statistics in 
the peripheral data at high pT is limited, the two most pe- 
ripheral bins are integrated for the fit in the high px range. 
The numerical values of T^f / as a function of centrality 
and of the different fit ranges are summarised in Tab. [3l 
An increase of the average Te/ / can be noticed when the 
fit is restricted to the low pT range, while a flatter trend 
is observed from the fit in the high px range. The differ- 
ence between the effective temperatures obtained in the 
two fit ranges is around 10 MeV (average over the three 
most central bins). This pattern indicates the presence of 
radial flow but it is rather small, as expected from the 
small coupling of the </> to the "pion wind" 5(),3T. 

Further quantitative insight has been gained studying 
the spectra in the framework of the blast wave model |52j : 



dN 



mTdrriT 



cx 



rdrniTlQ 



f Pt sinh p \ / tot cosh p 



To 



(2) 

where p = tanh Prif) and a linear velocity profile was 
assumed: (3r — Ps (f/R)- The slopes of the pt spectra de- 
fine lines in the Tq — Pt plane. Tq is varied in the window 
90-140 MeV. Transverse momentum distributions are then 
generated using eq. ([2]), with Pt values tuned in order to 
reproduce the observed Te// values. Fig. [T2l shows the re- 
sults for the (j) and also for the negative charged hadrons 
(essentially pions), rj, uj and p [51]. The study is performed 
for Pt > 0.4 GeV in order to exclude the rise in the low 
Pt region seen in the p excess spectra [S3]. In addition, 
to perform a comparison with meaningful statistical ac- 
curacy with the other particles, data were integrated in 
centrality, excluding the most peripheral events. A clear 
hierarchy in the freeze out parameters, as a consequence 
of the particle coupling to the medium is visible. The </> is 
the least coupled particle, while the p, being continuously 
produced in the medium via tt+tt" annihilation is maxi- 
mally coupled. The freeze-out component of the p receives 
the strongest radial flow boost, with a Te// > 300 MeV, 
50 MeV above the A l5Tll53l. 



4.4 Yield 

The average number of (p mesons per interaction (</>) is 
evaluated using the J/V' as a reference process. Once cor- 
rected for anomalous and nuclear absorption, {J/ip) scales 
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Fig. 12. Blast wave analysis results for 77, p, </> and nega- 
tive hadrons (pions). Fits performed in the range 0.4 < pr < 
1.8 GeV. 



with the number of binary collisions and can be written 
as 

(3) 



Here a 



j/ip,fj.fj. _ 

NN 



5.64 ± 0.1 nb is the J/ip cross section in 



nucleon-nucleon collisions multiplied for the J/^ — > pp 
branching ratio. It was derived by the systematic study of 
J/ij; production in pA collisions performed at 400 and 450 
GeV, rescaling to 158 GeV and extrapolating to full phase 
space [54j- The estimated nucleon-nucleon inelastic cross 
section at = 17.2 GeV is ctnn = 31.7 ± 0.5 mb [55] . 
((/)) can be evaluated from the experimental ratio: 



(0) 



N 



/( 



A 



X £ 



rec X J I 



(4) 



where N^^j/^^ is the number of (p {J/ip) mesons observed 
in a given centrality bin, A0(j/^,) x et'^^J^^ are the accep- 
tance and reconstruction efficiency, BRj^^ is the (j> pp 
branching ratio and ffj'^ incorporates the J/ip nuclear and 
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Fig. 13. Centrality dependence of the Te// parameter (full circles) compared to NA49 (open circles) and to NA50 (open 
squares). Left panel: fit performed in the NA49 range (0 < pt < 1-6 GeV). Right panel: fit performed in the NA50 range 
(PT > 1.1 GeV). 



anomalous absorption factors. The same centrality bins, 
selected using charged particle multiplicity as described 
in sec. 13.31 were used for both the (p and the J/ip. Invok- 
ing lepton universality, the value 2.97 x 10^* measured 
in the (/) — > ee is used for the branching ratio, since it 
is known with better precision [55j . The nuclear absorp- 
tion was accurately determined by several measurements 
in pA collisions at 400 and 450 GeV [SI] . The anomalous 
absorption was determined from the J/t/j analysis of the 
same In-In data sample [SB^. The overall correction fac- 
tor comes from the product of the nuclear and anomalous 
absorption. One has to notice that recent preliminary p- 
A results from NA60 showed the nuclear absorption at 
158 GeV is larger than previously assumed [S7|. However, 
since the anomalous absorption is defined relative to the 
nuclear one, the overall correction /jy*^ remains largely in- 
sensitive to relative changes of these two factors. The sys- 
tematic uncertainty arising from the anomalous and nu- 
clear absorption is 11% [56 . This is the dominating source 
of systematic uncertainty in the analysis of the yield. 

By multiplying the ratio in eq. ^ by eq.(I3]), one ob- 
tains (0) as a function of the number of binary collisions 
Ncoii. The values of NcoU and of Npart for each centrality 
bin are obtained by fitting the multiplicity distribution 
with a Glauber model, as described in sec. 13.31 and sum- 
marised in Tab. [TJ 

Processes having different scaling properties as a func- 
tion of centrality, such as the (f> and the J/tp, in general 
give rise to different multiplicity distributions and this 
could introduce a bias when using wide centrality bins. 
The analysis was repeated with a much finer centrality 
binning and this potential bias was found to be negligible. 



Table 4. (j> yield as a function of centrality. The bin by bin 
systematic error is reported in the table. The measurements 
are also affected by a 13% systematic uncertainty independent 
of centrality. 



Npart 




14 


0.172±0.007±0.021 


41 


0.690±0.025±0.074 


78 


1.48± 0.04 ±0.15 


132 


3.10± 0.08 ±0.2 


177 


5.0± 0.17 ±0.2 



The multiplicity integrated in centrality is (0) = 1.411± 
0.026 ±0.14. 

The robustness of the measurement is checked repeat- 
ing the analysis integrated in centrality, with a different 
method having largely independent systematics. The (p 
cross section is determined with the formula: 



f-trigf-DAQ BR^^^ 

where L is the integrated luminosity, e^v is the event se- 
lection efficiency, etrig is the trigger efficiency, euAQ is the 
efficiency related to the DAQ live-time. The integrated lu- 
minosity is given by L = NincNtgt, where Nine is the num- 
ber of incident ions as measured by the ZDC and Beam 
Tracker detectors, while Ntgt is the effective number of in- 
dium atoms in the target system [58] . The latter is affected 
by a systematic uncertainty of 10%, because the transverse 
profile of the beam is wider than the targets (except for 
the first one which had a radius of 0.6 cm). The trigger 
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Fig. 14. Te/ / for central collisions in several collision systems 
at 158 AGeV. Full circle: (f) ^ fifi in In-In (NA60); full upward 
triangle: K^K^ in p-p (NA49); full downward triangle: 
<j)^ K+K- in C-C and Si-Si (NA49); open circle: cj, ^ K+ K- 
in Pb-Pb (NA49); open square: — > in Pb-Pb (NA50); open 
triangle: (f> K+ K' in Pb-Pb (CERES). 



efficiency was determined with special runs taken with- 
out the dimuon trigger. It is 81% and is affected by a 6% 
systematic error. 

The average yield is then obtained dividing the cross 
section by the fraction of total In-In inelastic cross sec- 
tion, corresponding to the cut Nch > 4, as estimated 
with the Glauber Monte Carlo (see sec. 13. 3p . The total 
inelastic cross section calculated with the Glauber model 
is (jlnlr — 4-8 b, and has a 10% uncertainty, mostly due 
to the uncertainties on the nuclear density profiles. The 
centrality integrated (p cross section is cr^ = 6.83 ± 0.08 b. 
This leads to (</>) = 1.73 ± 0.02 ± 0.27, in agreement with 
the result based on the standard method, taking into ac- 
count only the unrelated systematic errors in the two ap- 
proaches. 

The results as a function of centrality, obtained by 
taking the average of the two methods weighted by the 
systematic error, are summarised in Tab. SI The ratio 
(0) /Npart is seen to increase with the collision central- 
ity indicating the presence of an enhancement as in other 
collision systems. 



4.5 Discussion of the results on T^jf and yield 

As mentioned in the introduction, the discrepancies in ab- 
solute yields and T^f f values observed by NA49 and NA50 
became historically known as the (p puzzle [I4lf22j . To ex- 
plain those differences, it was originally proposed that in- 
medium effects and kaon absorption or rescattering could 
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Npart 

Fig. 15. {(j>) /Npart in full phase space in several collision sys- 
tems at 158 AGeV. Full circles: (j) ^ fifi in In-In (NA60); full 
upward triangle: — » K" in p-p (NA49); full downward 
triangles: cf> — > K~ in C-C and Si-Si (NA49); open circles: 
(j) K^K~ in Pb-Pb (NA49); open squares: <;/>-» in Pb-Pb 
(NA50); open triangle: <t> K+ K' in Pb-Pb (CERES); open 
diamond: (/) ^ ee in Pb-Pb (CERES). Continuous line: NA50 
points extrapolated to full phase space with Te// as measured 
by NA49, dashed lines: uncertainty on these points. Square 
brackets: systematic error (the error bar of the Pb-Pb data in- 
cludes statistical and systematic error). C, Si and peripheral 
In points slightly displaced to improve readability. 



prevent the reconstruction of in-matter KK decays, 
in particular at low transverse momentum, while the 4> 
mesons decaying in the lepton channel would not be af- 
fected. This would lead to reduced yields and enhanced 
Teff values in the KK channel as compared to the /i/i 
channel. In this section, all the measurements performed 
at the SPS, including the new one in In-In, are compared 
and the hypothesis of the physical mechanism described 
above is critically reviewed, considering the facts in favour 
of and against it. 

Fig. [13] shows the NA60 T^f / obtained performing the 
fits in the pr windows corresponding to the coverage of 
NA49 (left) and NA50 (right). As mentioned in sec. [43l 
the difference of the temperature in the two pt windows in 
In-In is around 10 MeV, indicating a modest radial flow 
effect. Although a similar pattern is observed in Pb-Pb 
collisions - a stronger increase at low px (NA49) and a 
flattening at high pT (NA50) - the maximum difference 
is, for the most central Pb-Pb collisions, about 70 MeV. 
This difference seems not to be consistent with the Tf.ff 
variation implied by radial flow as observed in the NA60 
data. 

In the common fit range, the NA60 T^f / values in In-In 
are ~ 10 MeV larger than those measured by NA50, while 
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there is a seeming agreement between NA60 and NA49 in 
the common fit range. However, the NA49 errors are large 
in this region and, in addition, a direct comparison of the 
values of T^f / between NA60 and the Pb-Pb experiments 
for the same Npart is not straightforward, due to the dif- 
ferent geometry of the coUision (the relative weight of the 
thin nuclear halo is higher in peripheral Pb-Pb collision 
than in a smaller system). 

The temperature variation in different collision sys- 
tems, when restricting to central collisions only, is shown 
in Fig. [m This picture includes also the NA49 results 
for C-C, Si-Si and the CERES result in the kaon chan- 
nel for central Pb-Pb (the dielectron measurement has a 
very large error). The NA49 and CERES Pb-Pb points are 
above the In-In point by > 30 MeV. Part of this difference 
can be related to radial flow, since in Pb-Pb a larger life- 
time of the fireball would lead to a larger effect in central 
collisions. However, the T^ff measured for central Pb-Pb 
collisions is ^ 290 MeV, very close to that of the pro- 
ton [59]. The blast-wave analysis of the NA60 data, on 
the other hand, shows a clear hierarchy in the freeze-out, 
with the p having a maximal flow due to its coupling to 
pions and the being the least coupled meson fFig. [T^ . 
The temperature of the freeze-out component of the p in 
In-In reaches 300 MeV or more, 50 MeV above the 0. One 
would therefore expect an even larger difference in Pb- 
Pb collisions between the T^f / of the cj) and of the proton, 
which is strongly coupled to the pions via the A resonance 
and with a mass larger than the p. This suggests that the 
high Tef f seen by NA49 and CERES in the kaon channel 
cannot be ascribed entirely to radial flow. 

The (j) enhancement in In-In is compared to previous 
SPS measurements [THl[Tni[lIl[lS] in Fig. M 

An unambiguous comparison to the NA50 result in full 
phase space is not possible, due to lack of consensus on 
the value of the Te// parameter in Pb-Pb collisions, but 
in general the NA50 enhancement is considerably higher 
than what is seen in In-In at any centrality. An extrapola- 
tion to fuU pT with the NA50 T^jf values (T ~ 220 - 230 
MeV independent of centrality) leads to values higher by 
a factor 2 or more with respect to In-In. Even extrapo- 
lating with the NA49 T^jf values (T ~ 250 - 300 MeV) 
would lead to values significantly higher (as shown by the 
continuous line in Fig. [151 the dashed lines represent the 
propagated uncertainty on these values). 

The NA49 measurement^ show that the {(p) breaks 
Npart scaling and that (</>) /Npart, for central collisions as 
a function of the colliding system, already saturates for 
Si-Si. The enhancement in In-In is closer to NA49 than 
NA50, but it is not in agreement either. The two In-In 
most peripheral points are in agreement with C-C and 
Si-Si. However, for increasing centrality there is no sat- 
uration in In-In, while the NA49 Pb-Pb points show a 
flatter trend, which remains consistently smaller at any 

^ NA49 has characterized centrality of the collision using 2 
quantities: the "number of wounded nucleons" and the "num- 
ber of participants" . The one directly comparable to our A'part 
is the number of wounded nucleons. The values for the different 
Pb-Pb bins are taken from Ref. [59j . 



centrality. One could ask whether the increase of the yield 
ratio NA60/NA49 as a function of centrality could be as- 
cribed to a suppression mechanism present in Pb-Pb that 
tends to decrease the overall enhancement. Restricting the 
comparison among In-In and Pb-Pb central collisions, the 
dimuon NA60 measurement remains larger than the kaon 
Pb-Pb points. 

Summarizing, the large NA49/NA50 difference in cen- 
tral collisions for T^/f (70 MeV) and the yield (a factor 4), 
was hardly accounted for by theoretical models. A differ- 
ence, though smaller, is also present in central collisions 
in Tgff and yield, when considering the NA60 muon mea- 
surement. This difference remains difficult to explain and 
hints for the presence of a physical effect. Whether this 
is present also in In-In, or peculiar only to Pb-Pb, would 
require the comparison to ^ KK in In-In. Preliminary 
NA60 results of K^K~ suggest that there are no 
differences within the errors among the two channels in 
In-In [60 . 

One has to notice that the CERES dielectron and kaon 
points for central Pb-Pb are in agreement and in addition, 
both the NA50 Teff and yield do not agree with NA60 
either. Thus, at the moment it is not possible to include 
into a coherent picture all the measurements performed so 
far in the different collision systems. 

4.6 Mass and width 

As mentioned in the introduction, the spectral function 
of the (f> could be modified in the medium, resulting in 
changes of its mass and width. The lifetime of the (p is 
44 fm, larger than the estimated fireball lifetime (~ 7 fm). 
Any medium modified component would therefore sit un- 
der a large unmodified peak produced after freeze-out. In 
order to distinguish a possible in-medium component from 
instrumental effects, a careful tuning of the Monte Carlo 
simulation is required. This was not strictly necessary in 
previous analysis of the low mass region excess [ST I IHT j , as 
this structure, mostly due to tttt p* —f up, represents 
a large fraction of the total yield (the lifetime of the p is 
- 1.3 fm). 

Uncertainties in the detectors alignment, materials, en- 
ergy loss compensation and precision of the magnetic field 
maps lead to residual differences between the Monte Carlo 
and the real data. These were compensated by fine-tuning 
the Monte Carlo against the peripheral data, where no in- 
medium effects are expected. The tuning was performed 
considering the whole mass spectrum below 1 GeV and 
required a global 0.3% ± 0.1% mass shift and a 4% ± 1% 
change in the mass resolution (previous NA60 analyses 
only included the mass shift correction). The uncertain- 
ties in this tuning are the main source of systematic error 
for the results discussed in this section. 

In order to look for an in-medium component around 
the peak, the excess was isolated using the same sub- 
traction procedure already extensively used for the low 
mass region (61]. The w, and r] Dalitz were subtracted 
from the total data using the Monte Carlo shapes and fix- 
ing the yield using solely local criteria; the uj Dalitz and 
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Fig. 16. Excess mass spectrum for data integrated in centrality, obtained without (left) and with (right) the Monte Carlo fine 
tuning. Red circles: data; black triangles: excess data; lines: individual processes. 



the T] two body were then subtracted by binding them 
to these processes via branching ratios. The quality of the 
subtraction can be judged in Fig. [TBI which shows, for data 
integrated over centrality, the excess before and after the 
fine tuning. While the tuning is essentially constrained by 
the narrow cj) peak, it improves the subtraction over the 
whole mass spectrum: the evident artifacts visible in the 
left panel around all the narrow resonances (ry, w and (j)) 
disappear after the tuning. The continuum left after the 
subtraction under the (j) peak is smooth: no evidence for 
a medium modified can be seen in this figure. Indeed, 
an in-medium broadened component would manifest itself 
as broader tails extending to the sides. In this case, the 
subtraction of the vacuum component should not leave a 
smooth continuum, but rather a broad structure or a dip, 
which is not observed. 

As mentioned in the introduction, the KEK-PS E325 
experiment reported evidence for a large excess to the left 
of peak in p-Cu at 12 GeV for (3^ < 1.25, corresponding 
to a produced (j) moving very slowly inside the heavy target 
nucleus at rest. This situation is not directly comparable 
to the high energy regime discussed in this paper. Here, a 
fireball is produced with particle momenta boosted in the 
forward direction, therefore the relevant variable to select 
(f) mesons mostly sensing the medium is transverse mo- 
mentum. Moreover, in a high energy heavy ion collision, 
as opposed to p-A, the can be produced not only at 
chemical freeze-out, but also in the hot hadronic medium 
via kaon reactions. Fig. [T7] shows the region around the 
(j) mass peak for different px and centrality selections: the 
left panel shows the two most central bins and for pT < 0.6 
GeV, while the right panel shows the data integrated in 
centrality for px < 0.2 GeV. The continuum from the 



subtraction is smooth also in this case. Furthermore, a fit 
was performed using the Monte Carlo shapes for the dif- 
ferent processes. The continuum was described with the 
model of Ref.|llj. which includes, besides the in- medium 
p, multi-pion and partonic processes relevant in this mass 
region. The model of Ref . [62] , having different weights for 
the multi-pion and partonic processes, does not change 
the conclusions, since it produces a similar yield below 
the (j) peak. The /ndf is ~ 1 in all cases. There is no 
evidence of any mass shift or broadening, even at high 
centrality and low pT , where one would expect to observe 
the strongest in-medium effects. 

In order to estimate quantitatively the pole mass and 
total width of the (p peak, data were fitted with the con- 
volution of a relativistic Breit-Wigner and of a function 
parametrizing the detector response, on top of a linear 
background. In a first step, the parameters of the detector 
response were adjusted on the Monte Carlo simulation. 
The data are then fitted leaving the mass and width of 
the Breit-Wigner as free parameters. The quality of the 
fit is very good, as shown in the left panel of Fig. [T8l 
corresponding to the semi-central bin, Npart = 133. The 
results the mass and width as a function of centrality are 
shown in the right panel of Fig. [181 The systematic er- 
ror was estimated changing the fit range, the background 
function and propagating the uncertainties in the Monte 
Carlo tuning. The mass and width are compatible with 
the PDG values and independent of centrality, although 
our sensitivity on the width is limited by the mass resolu- 
tion, which is much larger than the natural width (this is 
reflected in the large statistical error in the right panel of 
Fig. [HI). 
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Fig. 17. <j> mass region in the two most central bins for pr < 0.6 GeV (left) and all centralities, pr < 0.2 GeV (right). 
Red circles: data; black triangles: excess data; lines: individual processes (RH: model of Ref. [H]); black crosses: result of fit 
performed using the Monte Carlo shapes. 




Fig. 18. Left panel: fit of the (j) mass peak, semi-central bin {Npart = 133) with a relativistic Breit-Wigner convoluted with 
the detector response, on top of a linear background. Right panel: mass (top) and width (bottom) of the as a function of 
centrality. Shaded box: systematic uncertainty common to all centrality bins. 
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While the data can be described by a vacuum 0, we 
set an upper limit on a modified component using the 
only existing prediction for in-medium production of <j) 
in 158 AGeV In-In collisions [ll]. This corresponds to 
an in-medium spectral function folded over the phase- 
space fireball evolution up to freeze-out and appears sig- 
nificantly broadened (> 80 MeV), so that the data should 
have sufficient sensitivity - in terms of mass resolution 
- to discriminate it. The region around the (p peak was 
fitted adding an in-medium cf) component to the unmod- 
ified one. In the centrality integrated data, the fraction 
in-medium/ vacuum 4> is compatible with zero within one 
sigma with an upper limit of 0.08 ± 0.03 (syst) at 95% 
confidence level. In the two most central bins the fraction 
remains compatible with zero within one sigma and the 
upper hmit increases to 0.15 ± 0.03 (syst), due the larger 
statistical uncertainty. The fraction predicted by Ref. [11] 
for semi-central data (roughly equivalent to centrality in- 
tegrated data) is ~ 4%. The limitation in statistics pre- 
vents us from setting a meaningful upper limit at high 
centrality and low pT- 

5 Conclusions and outlook 

In this paper a complete characterization of — > /i/z pro- 
duction in In-In collisions has been provided by measuring 
the differential spectra, yield, mass and width. 

The rapidity distribution width does not depend on 
centrality within errors, and is found to be compatible 
with previous SPS measurements. The decay angle distri- 
butions show no evidence for polarization in any of the 
reference frames used, independent of centrality and pr- 
This confirms the lack of anisotropy previously reported 
by NA60. This is a necessary though not sufficient condi- 
tion for particle production from a thermalizcd medium. 
Indeed, the centrality dependence suggest that the po- 
larization is zero or very close to zero also in peripheral 
collisions, close to the limit of p-p like interactions. 

The study of the transverse momentum distributions 
shows a small dependence of the Te// parameter on the 
Pt window considered, to be ascribed presumably to ra- 
dial flow. A blast wave analysis shows a clear hierarchy of 
the freeze-out, with the </) being the least coupled particle 
to the medium. The difference between Tg// measured in 
Pb-Pb collisions by NA49 and NA50 is much larger and 
cannot be explained by the different pt windows explored. 
The Tf,f / value for most central collisions in In-In remains 
significantly lower than the kaon point in Pb-Pb measured 
by NA49 and CERES. 

The yield measured by NA50 at high px, when ex- 
trapolated to full Pt, exceeds the NA60 data points for 
the same number of participants and seems to lie out- 
side of the general systematics. The NA60 yield, on the 
other hand, remains consistently larger than the NA49 
and CERES Pb-Pb measurement. 

The new measurements in T^ff and yield, when com- 
pared to the kaon channel in Pb-Pb, suggest that there 
may be a physical effect leading to a difference in the two 
channels. Further clarification may come from the NA60 



measurement of the (p —t K~ channel in In-In colli- 
sions, which is currently in progress and will be reported 
in a forthcoming paper. At the moment, it is not possible 
to include all the existing SPS measurements in a coherent 
framework. 

The mass and width were studied as a function of 
centrality and transverse momentum. No evidence for in- 
medium modifications was found. Our data were also com- 
pared to the only existing prediction for in-medium cj), pro- 
duced in In-In collisions at 158 AGeV. The yield of such 
a component is compatible with zero, with an upper limit 
of about 8% for the centrality integrated data. 
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